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ABSTRACT 

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds 
Dark Nebulae with a resolution of w 2' are reported. These sources are selected from the 
SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched 
to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with 
radio observations at lower frequencies with matched Mv-plane coverage is made, in order to 
search for any anomalous excess emission which can be attributed to spinning dust. Possible 
emission from spinning dust is identified as a source within a 2 arcmin radius of the Scuba 
position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency ra- 
dio emission. We find five sources which show a possible spinning dust component in their 
spectra. These sources have rising spectral indices in the frequency range 14.2-17.9 GHz with 
a\ll = -0.7 ±0.7 to -2.9 ±0.4, where 5 oc v"". Of these five one has already been reported, 
LI 1 1 1, we report one new definite detection, L675 (16 a), and three new probable detections 
(L944, LI 103 and L1246). The relative certainty of these detections is assessed on the basis 
of three criteria: the extent of the emission, the coincidence of the emission with the Scuba 
position and the likelihood of alternative explanations for the excess. Extended microwave 
emission makes the likelihood of the anomalous emission arising as a consequence of a radio 
counterpart to a protostar or a proto-planetary disk unlikely. We use a 2 arcmin radius in order 
to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third 
criterion is used in the case of LI 103 where a high flux density at 850 /xm relative to the FIR 
data suggests a more complicated emission spectrum. 

Key words: Radiation mechanisms:general - lSM:general - ISMxlouds - Radio contin- 
uum:general 



1 INTRODUCTION 

Modem theoretical models of the infrared emission from dust in 
the diffuse ISM require a large proportion of starlight energy to be 
reradiated by ultrasmall grains (~20%; Draine 2003). Comparison 
of the observed spectra with theory suggests that a large fraction 
of these small grains must be in the form of Polycyclic Aromatic 
Hydrocarbons (PAHs). It is supposed that in addition to their vibra- 

* We request that any reference to this paper cites "AMI Consortium: Scaife 
et al. 2009" 

t E-mail: as595@mrao.cam.ac.uk 



tional emission, these small grains will also rotate (Erickson 1957; 
Ferrara & Detmar 1994; Draine & Lazarian 1998). The emission 
from this rotation is found, using simple thermodynamical argu- 
ments, to be in the cm-wave regime of the electromagnetic spec- 
trum. The details of the rotational excitation and damping have 
been examined in detail by Draine & Lazarian (1998) and their 
models have been found to be in good argreement with observa- 
tional results. 

Although the anomalous emission was originally seen as a 
large-scale phenomenon (Leitch et al. 1997, see also Kogut et al. 
1996) in CMB observations it has also been seen to occur in a num- 
ber of different types of astronomical objects, such as Hll regions. 
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dark clouds and photodissociation regions (Casassus et al. 2006; 
Dickinson et al. 2007; Scaife et al. 2007; AMI Consortium: Scaife 
et al. 2009). It is often found to be correlated with thermal dust 
emission, but this is not always the case (Casassus et al. 2008). A 
growing number of observations have detected dust-correlated cm- 
wave emission which is inconsistent with free-free emission. How- 
ever, these few observations are not yet sufflcient to make any def- 
inite observational statement about the enviromental factors which 
influence whether emission from spinning grains will be seen to- 
wards a certain object. We hope that the observations presented 
here will contribute significantly to improving this knowledge. 

In this paper we have imaged fourteen dark molecular clouds. 
Our sample is selected from the submillimetre survey of dark 
molecular clouds performed using the Submillimetre Common- 
User Bolometre Array (SCUBA) on the James Clerk Maxwell Tele- 
scope (JCMT) by Visser, Richer & Chandler (2001). Of these four- 
teen only three are currently known to contain protostars (L944, 
L1014 and L1246; Visser et al. 2001; Young et al. 2004), with the 
remainder thought to be examples of pre-stellar cores. 

The dark molecular clouds in this sample were taken from the 
Lynds dark nebula catalogue (Lynds 1962) of 1802 optically se- 
lected dark clouds from the POSS plates. The SCUBA jiggle map- 
ping function used to observe the sample restricted the field of view 
to 2.3 arcmin diameter This constraint led to a selection criterion of 
a maximum semi-minor axis of < 1.1 arcmin for each cloud, which 
was determined from their optical size as listed in Parker (1988). 

The Lynds catalogue is divided into opacity classes from 1 to 6 
of which 6 is the most opaque and includes 147 clouds. The clouds 
observed in the SCUBA survey are all class 6 (Ay ^ 5 mag; Lynds 
1962) and should therefore provide the highest column density of 
dust when searching for anomalous emission from spinning grains. 

We have observed all clouds in the sample of Visser et al. 
(2001) which fall in the available declination range of the Ar- 
cminute Microkelvin Imager (AMI) Small Array (SA) instrvunent 
in Cambridge. The small angular size of these clouds makes them 
well matched to the size of the AMI SA synthesized beam. We iden- 
tify excess microwave emission from the SCUBA clouds as any 
source which falls within a 2 arcmin radius of the SCUBA position 
and which shows an excess of emission relative to lower frequency 
radio data. In addition, any source which satisfies these two criteria 
but has a falling spectrum across the AMI band will be rejected as a 
spinning dust candidate. This allows us to distinguish those sources 
which may be optically thick at lower radio frequencies, but have 
spectra which turn over before the AMI band. 

The paper is structured as follows: Section 2 describes the tele- 
scope and the observations. Section 3 details the methods wc use to 
compare our data to lower frequency observations, Section 4 dis- 
cusses each object individually and Section 5 discusses the sample 
as a whole and states our conclusions. 



2 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(Small Array) and one of eight 13 m antennas (Large Array), both 
sited at Lord's Bridge, Cambridge (AMI Consortium: Zwart et al. 
2008). The telescope observes in the band 13.5-17.9 GHz with 
cryostatically cooled NRAO indium-phosphide front-end ampli- 
fiers. The overall system temperature is approximately 25 K. The 
radio frequency is mixed with a 24 GHz local oscillator, down con- 
verting to an intermediate frequency band of 6-10.5 GHz. Amplifi- 
cation, equalization, path compensation and automatic gain con- 



trol are then applied to the IF signal. The correlator is an ana- 
logue Fourier transform spectrometer with 16 correlations formed 
for each baseline at path delays spaced by 26 rrmi. In addition both 
real and imaginary correlations are formed by use of 0° and 180° 
hybrids. From these, eight 0.75 GHz bandwidth channels are syn- 
thesized. In practice, the two lowest frequency channels (1 & 2) are 
not generally used due to a low response in this frequency range 
and interference from geostationary satelUtes. 

Observations of fourteen dark molecular clouds from the 
Lynds dark nebula catalogue were made with the AMI SA tele- 
scope during the period 2007 January - May. The AMI LA was not 
used in this work. Here we present data for each of these regions, 
which are listed in Table 1 . These targets were chosen from the sub- 
millimetre survey of Visser et al. (2001) as small molecular clouds, 
well matched to the AMI synthesized beam, in a declination range 
consistent with the capabilities of the AMI at a latitude of 52° N. 
Each pointing was observed at least twice for at least eight hours, 
with the exception of LI 185 which was observed only once due to 
poor weather. 

Data reduction was performed using the local software tool 
REDUCE. This applies both automatic and manual flags for inter- 
ference and shadowing and hardware errors. It also apphes phase 
and amplitude calibrations; it then Fourier transforms the correla- 
tor data to synthesize frequency channels before output to disk in 
uv FITS format suitable for imaging in AIPS. 

Flux cahbration was performed using short observations of 
3C286 near the beginning and end of each run. We assumed I-^Q 
flux densities for this source in the AMI SA channels consis- 
tent with the frequency dependent model of Baars et al. (1977), 
~ 3.3 Jy at 16 GHz. As Baars et al. measure I and AMI SA mea- 
sures I-i-Q, these flux densities include corrections for the polari- 
sation of the calibrator source derived by interpolating from VLA 
5, 8 and 22 GHz observations'. A correction is also made for the 
changing intervening air mass over the observation. From other 
measurements, we find the flux calibration is accurate to better than 
5 per cent (AMI Consortium: Scaife et al. 2008; AMI Consortium: 
Hurley-Walker et al. 2009). 

The phase was calibrated using hourly interleaved observa- 
tions of calibrators selected from the Jodrell Bank VLA Survey 
(JVAS; Patnaik et al. 1992). After calibration, the phase is gener- 
ally stable to 5° for channels 4-7, and 10° for channels 3 and 8. The 
FWHM of the primary beam of the AMI SA is ^ 20' at 16 GHz. 
Channels 1 and 2 are generally discarded due to heavy satellite in- 
terference. 

Reduced data were imaged using the AIPS data package. Both 
the MEM and CLEAN deconvolution methods were employed, as 
outlined in the next section. CLEAN deconvolution was performed 
using the task imagr which applies a differential primary beam 
correction to the individual frequency channels to produce the com- 
bined frequency image. CLEAN deconvolution maps were made 
from both the combined channel set and for individual channels. 
The full set of CLEANed AMI maps is shown in Appendix A. With 
the exception of Fig. 2 all of the images presented in this paper 
were made using clean deconvolution and are not primary beam 
corrected. MEM deconvolution is also used on occasion in this pa- 
per, for reasons described in Section 3, and is implemented through 
the AIPS task vtess. The broad spectral coverage of AMI allows 
a representation of the spectrum between 14.2 and 17.9 GHz to be 
made independently of other telescopes and in what follows we use 
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Table 1. AMI Lynds cloud sample. Column [1] gives the source name, Column [2] the Right Ascension, Column [3] the declination, Column [4] Usts the r.m.s. 
specific intensity outside the primary beam for each map. Column [5] lists the optical Semi-major Axis of the cloud, Colunm [6] Usts the optical Semi-minor 
Axis of the cloud, Colunm [7] Usts the major axis of the synthesized beam. Column [8] lists the minor axis of the synthesized beam and Colunm [9] Usts the 
complex associated with each dark cloud. 



Name 


ra' 


Dec' 


Olms 


Semi-major Axis^ 


Semi-minor Axis- 






Associated 




(J20()0) 


(J20()0) 


^ beam ' 


(arcmin) 


(arcmin) 


(arcmin) 


(arcmin) 


complex 


L675 


19 23 52.6 


1 1 07 39 


0.17 


2.2 


0.6 


2.6 


2.4 




L709 


19 13 54.1 


16 26 58 


0.31 


1.1 


1.1 


2.5 


2.4 




L860 


20 03 18.1 


36 01 45 


0.21 


1.1 


1.1 


2.5 


2.1 


Cloud B 


L917 


20 39 55.1 


44 08 41 


0.27 


1.7 


0.6 


2.8 


2.1 


Cloud B 


L944 


21 17 40.8 


43 18 08 


0.27 


1.7 


1.1 


2.5 


2.3 


Cygnus Rift 


L951 


21 20 13.0 


43 31 30 


0.25 


2.2 


0.6 


2.5 


2.2 


Cygnus Rift 


L953 


21 21 25.5 


43 21 04 


0.55 


1.1 


0.6 


2.5 


2.4 


Cygnus Rift 


L1014 


21 24 01.7 


49 59 26 


0.38 


1.7 


0.6 


2.5 


2.1 


Cygnus Rift 


L1021 


21 21 39.3 


50 57 20 


0.40 


1.1 


0.6 


2.4 


2.1 


Cygnus Rift 


Ln03 


21 42 10.2 


56 43 44 


0.48 


2.8 


1.1 


2.4 


2.1 




Lllll 


2140 27.1 


57 48 10 


0.29 


1.7 


1.1 


2.4 


2.1 




L1166 


22 05 26.2 


59 33 38 


0.56 


1.1 


0.6 


2.6 


2.2 




L1185 


22 29 24.5 


59 09 23 


1.12 


2.8 


0.6 


2.8 


2.2 




L1246 


23 25 12.9 


63 39 30 


0.36 


2.8 


1.1 


2.5 


2.2 


Cepheus OB3 
complex 



' positions adapted from Visser et al. (2001) 
2 Parker (1988) 



the convention: 5 o= V " , where S is flux density, V is frequency and 
a is the spectral index. 



3 MULTI-FREQUENCY ANALYSIS 

The synthesized beam of AMI SA has FWHM rj 2-3 arcminutes. 
This makes comparison with low frequency surveys difficult as the 
GB6 survey at 4.85 GHz (Gregory et al. 1996) and the Effelsberg 
telescope surveys at 2.7 GHz (Fiirst et al. 1990) and 1.4 GHz (Reich 
et al. 1997) have larger beams. Moreover, flux losses are important 
in the context of the observations presented here as all of the sam- 
ple of dark clouds lie within the Galactic plane, where the emission 
on different scales is complex, and will be represented differently 
by telescopes with different effective wv-coverage. To constrain the 
spectral behaviour of these sources we employ two uv matching 
strategies which are outlined below. Similarly the information pro- 
vided on these regions at 100 flm by the IRAS satellite and its re- 
cently re-processed versions (Miville-Deschenes & Lagache 2005; 
Schlegel, Finkbeiner & Davis 1998) are at a resolution of ~ 4 ar- 
cmin, larger than that of AMI. Consequently we rely on the finer 
resolution results of the SCUBA observations to constrain the dust 
properties of the clouds. 

In order to obtain spectra for these sources, direct comparison 
with lower frequency surveys is not possible. We must therefore 
perform some operation on the data in order to match the observed 
angular scales. In the case that fcsD bo, where ^sd is the beam 
of the single dish telescope and /^o is the synthesized beam of the 
interferometer, the two sets of visibilities and their resulting maps 
can be compared directly with little loss of information, since the 
Fourier transform of the single dish beam, 6sd(M;V) 1 over the 
(«, v) range of interest. However, this is not true in the case that bsD 
approaches or is greater than hf). 

For nine of our fourteen Lynds clouds lower frequency radio 
observations are available at 1.4 GHz from the Canadian Galactic 
Plane Survey (Taylor et al. 2003; CGPS hereinafter). The CGPS is 



a combination of fully sampled interferometric data and single dish 
data, which provides a total power measurement of the sky with a 
beam of 60" -C ^ami — 150" (depending on declination). For the 
nine clouds covered by the CGPS we can make the approxima- 
tion that ^sD (M; v) 1 over the uv range of interest and compare 
«v-sampled CGPS data directly to our measured visibilities from 
AMI. However, in the case of the remaining five clouds we are lim- 
ited by the resolution of the available data and can only make a 
map plane comparison, after convolving the clean AMI map to a 
common resolution. This leads to two distinct analysis paths. 

The details of these data reduction methods are of importance 
when comparing flux densities from different instruments. This is 
especially the case when using interferometers. We include them 
here to emphasise the point that when we compare AMI observa- 
tions to those at lower radio frequencies we are comparing matched 
angular scales, unless stated otherwise. 

3.1 Image reconstruction 

To simulate the visibilities we take the single dish map and multiply 
by the primary beam of AMI, we then FFT the modulated map 
to obtain the data in the form of visibilities. These visibilities are 
sampled according to the measured visibilities in the matching AMI 
observation. The sets of sampled visibilities are each imaged and 
deconvolved separately to obtain our final simulated maps. 
Case 1: CGPS data, bso <: bo: We use data from the CGPS 
archive at 1.4 GHz for simulations. These data have a resolution 
of ^SD.FWHM ~ 60". We can examine the correlation between the 
two data sets in both the map plane and the uv plane. 
Case 2: non-CGPS data, Z?sd ^ ^o: For example, data from the Ef- 
felsberg 100 m telescope survey at 2.7 GHz. These data have a reso- 
lution of fcsD.FWHM = 4^3. In this case we cannot directly compare 
the visibilities in the uv plane, since the data is multiplied by the 
Fourier transform of the convolving beam, b^u- In the map plane, 
however, we can recover usable datasets for comparison when the 
MEM method of deconvolution is applied. 
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In brief, when deconvolving interferometric data using the 
CLEAN algorithm a recovered sky will have the form: 

^(^wWovered = [I{l,m) X B{l,m)]* bo{l,m) (1) 

where 1(1, m) is the true sky distribution, B{l,in) is the primary 
beam response and b(,{l,m) is the cleEin beam. In the case where 
we sample single dish data our input intensity distribution is no 
longer the true sky but instead a model sky, which is the convolution 
of the true sky with the single dish beam, l'{l,m) = I{l,m) *i'sD- 
Since this convolution is non-associative with the multiplication by 
the primary beam response, the two maps cannot be quantitatively 
compared. 

However, if instead we choose to deconvolve using an MEM 
method, such as the AIPS task VTESS, the primary beam response 
is accounted for in the fitting. MEM lets us recover: 

I {I, 'M)recovered = l' {l,m) *bo{l,m) 

= I{l,m)*bsu*bo{l,m). (2) 

Since convolution is associative with itself we can then convolve 
our own AMI data to match our sampled data. We do this by de- 
convolving the AMI data using VTESS to recover 1(1, m) *h[)(l,m) 
and convolving this recovered map with the single dish beam bgD- 
This will result in our comparison maps, both the sampled data 
and the convolved AMI data, having a clean beam of b(l,m) = 
\/bo(l,m)'^ + bsY)(l,m)'^. We note that this process will degrade the 
original AMI data, however it results in maps which have identical 
resolution and contain the same angular scales. 

3.2 Flux extraction 

The morphology of the sources in these observations is often not 
well described by a Gaussian fit. We therefore estimate their flux 
densities by removing a tilted plane fitted to the local background 
and integrating the remaining flux, see for example Green (2007). 
We do this by drawing a polygon around each source and fitting a 
tilted plane to the pixels around the edge of the polygon. Where an 
edge of the polygon crosses a region confused by another source 
the background is subjective and we omit this edge from the fit- 
ting. Example polygons are shown in e.g. Fig. 2 with omitted edges 
shown as dashed lines. Where we believe there may be some con- 
fusion as to the aperture selected in what follows we have depicted 
an example aperture in the same way, see the figure captions for 
details. Since this method is dependent to some degree on the aper- 
ture selected around the source we perform five fits for each object 
changing the aperture slightly each time, and take the flux density 
as the average of these fits. The vertices of these apertures are listed 
in Appendix A2. We denote the variance of these fits as cr^j and in- 
clude it in the final error on the calculated flux density as: 

a = ^(t2„, + (0.055)2 + (3) 

Here we represent the r.m.s. fluctuations outside the primary beam 
measured from the map as arms and the flux density of the source 
as 5. This calculation assumes a conservative 5 percent error on the 
flux calibration. 

We note that images made using CLEAN are often system- 
atically different to those made using MEM through VTESS. The 
major, although not only, difference arises due to the positive def- 
inite nature of the model assumed in VTESS. This has the effect 
of raising the background level in images reconstructed by VTESS 
relative to those made using CLEAN. Consequently it is not rec- 
ommended to use these MEM reconstructions for comparison with 



images restored with CLEAN. Our method of flux extraction acts 
to mitigate this effect through subtracting a tilted plane. However, 
we note that we have not used MEM images for quantitative com- 
parison except in the case of L675 where we have compared like 
with like. In later sections we have sometimes included a flux den- 
sity recovered using an MEM map to enhance a flux spectrum, but 
this is purely for illustrative purposes and will be expUcitly com- 
mented upon in the text. We reiterate: Case (1) employs CLEAN 
deconvolution. Case (2) employs MEM deconvolution. 

3.3 Sub-mm data 

The SCUBA data of Visser et al. (2001;2002) in conjunction with 
IRAS data available from the literature allows us to place con- 
straints on the thermal dust spectrum of these objects and their dust 
temperature, T^. 

Following Andre, Ward-Thompson & Barsony (1993) we use 
a modified greybody spectrum of the form 

Sy=By{Td){\-er'^)as,y (4) 

to fit these data, where Bv(Tj) is the Planck spectrum for a tem- 
perature at a frequency v. Ty is the optical depth of the cloud, 
generaUy assumed to be proportional to v^-^, and £2^ y is the extent 
of the source at a given frequency. The inclusion of tls,v allows us 
to account for the different source sizes seen by IRAS and SCUBA. 



4 roENTIFICATION OF RADIO COUNTERPARTS 

In this section we describe the radio emission seen towards the 
Lynds dark nebulae listed in Table 1 between 14.2 and 17.9 GHz. 
Sources are listed in Table 2. In what follows we identify potential 
excess microwave emission from the SCUBA clouds as any source 
which falls within a 2 arcmin radius of the SCUBA position and 
which shows an excess of emission relative to lower frequency ra- 
dio data. In addition, any source which satisfies these two criteria 
but has a falling spectrum across the AMI band will be rejected as a 
spinning dust candidate. This allows us to distinguish those sources 
which may be optically thick at lower radio frequencies, but have 
spectra which turn over before the AMI band. All errors are quoted 
to 1 a. 

4.1 L675: 

L675 shows a large radio counterpart at 16 GHz. This region lies 
outside the the limits of the CGPS data, however it is covered by 
NVSS at 1.4 GHz, Effelsberg at 2.7 GHz and the GB6 survey at 
4.85 GHz. These data show no counterpart to the object seen in the 
higher frequency AMI map, see Fig. 1. 

Since there is no coverage of this region by the CGPS we use 
our analysis Case (2) to determine the spectrum of this object. Sam- 
pling data at 2.7 GHz and 4.85 GHz from the Effelsberg 100 m tele- 
scope' and the GB6 survey^ respectively, and correcting for the 
beams as described in Section 3 it is immediately noticeable that 
the morphology of this object is very different at 16 GHz to that of 
the lower frequencies: see Fig. 2. A correlation analysis of the maps 
shows that the AMI and Effelsberg data have a Pearson coefficient 
of only 0.25, and AMI and GB6 have a slightly higher correlation 

' http : //www . mpif r-bonn . mpg . de/old_mpifr/ survey . html 
^ http : / /skyview . gsf c . nasa . gov/cgi-bin/query . pi 



AMI observations of Dark Nebulae 5 



a ii°05'- 




19°24'°30' 24'°00= ZS^ZO' SS^OO" 

Right Ascension (J2000) 
(a) 



^ iino'- 



a ir-os'- 




19°24"30= 24='00= 23"30' 
Riglit Ascension (J2000) 
(b) 



Figure 2. L675: (a) AMI 16 GHz combined channel data convolved to mtach the resolution of Effelsberg 2.7 GHz and (b) Effelsberg 2.7 GHz data sampled to 
match the uv coverage of AMI. Both maps have a resolution of 5'A x 5(2, shown as a fiUed ellipse in the bottom left comer, and have been deconvolved using 
MEM. The AMI primary beam is shown as a circle and the pointing centre as a cross. These maps are primary beam corrected. Contoiffs are 10 percent from 
the 50 percent level. Example flux extraction apertures are shown where a dotted line indicates an edge not used for fitting a titled plane base-level, see text for 
details. 
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19''24='30' 24"00' 23"30= 23"00' 

Riglit Ascension (J2000) 

Figure 1. L675: CLEANed combined channel AMI data at 16 GHz are 
shown as contoiffs. Contoiffs are at 3, 6, 12, 18 ffnns etc. NVSS data are 
shown as greyscale. The AMI primary beam FWHM is shown as a soUd cir- 
cle and the pointing centre as a cross. The AMI synthesized beam, 2f6 x 2'A, 
is shown as a filled ellipse in the bottom left comer. This map is not primary 
beam corrected. 

coefficient of 0.38. The correlation of GB6 and Effelsberg is 0.56. 
However, the correlation between the AMI map and appropriately 
sampled IRIS 100 jUm data^ is 0.82, demonstrating that the 16 GHz 
emission follows the infra-red emission much more closely than 
the low frequency radio data. Since the object is not detected in 
the lower frequency maps we take the intensities at the position 
of the cloud at those wavelengths and regard them as upper limits 
on any emission. These values are 52.7 ^ 6.2 mjy beam"' for Ef- 
felsberg data and 54 gs < 3.8 mjy beam"' for the GB6 data. The 
spectral index, cHj 7^ ^ 0.84, is consistent with a region of non- 
thermal emission and predicts a flux density of 1.4 mjy at 16 GHz. 
At the same resolution, 5'.A x 5'.2, the AMI emission has a flux den- 
sity of 5i6 = 20.2 ± 1 .02 mJy. This indicates an excess of 16.4 mJy 



(16. 1 (7) at 16 GHz relative to a flat spectral index from 4.85 GHz or 
an excess of 18.8 mJy (18.4 <t) relative to the predicted flux density 
from Oj I^, assuming a source size equivalent to the beam at lower 
radio frequencies. Note that the noise level in the original GB6 data 
at this declination is ~ 5 mJy beam^' (Gregory et al. 1996) and in 
the original Effelsberg data is ^ 8mJybeam^' (Ftirst et al. 1990) 
placing all the structure within the AMI primary beam below the 
2a noise level at both frequencies. 

Two further sources are evident within the primary beam of 
the AMI field. Of the two point sources. Source "A" lies almost ex- 
actly on the FWHM point of the beam to the north-west of L675. It 
appears in the six AMI channels and Effelsberg data with approx- 
imately the same flux density: 52,7 = 4.3 mJy. The best fit spec- 
tral index from these data is a = 0.00 ±0.15 . We identify this 
source as a region of optically thin thermal emission. Note that 
the derived flux density of this source at 2.7 GHz lies below the 
noise level of the original Effelsberg data and consequently we 
treat it as an upper limit. Source "B", see Fig. 2, appears both in 
the NVSS map as a compact object and in the Effelsberg data as a 
more extended source, with flux densities of 5i.4 = 40.0 ± 1.6 mJy 
and 52.7 = 13.1 ±4.9 mJy respectively. This indicates a steep spec- 
trum source and, although there may be a small amount of flux at 
16 GHz as indicated by the morphology of the main object, there is 
no identifiable discrete object at this frequency. 

Although there is no object associated with L675 in the IRAS 
point source catalogue (hereinafter IRAS PSC) it was part of a 
deep IRAS photometry catalogue (Clemens, Yuo & Heyer 1991; 
hereinafter CYH91). This catalogue yielded fluxes for L675 at all 
four IRAS frequencies and a size of Cls = 39 arcmin-^, which corre- 
sponds well to the AMI data (note that the source size from CYH91 
corresponds to all flux down to the fitted background, not a Gaus- 
sian FWHM). 

A greybody model, as described in Section 3.3, where Ty °= 
v''^, fitted to the 100/60 flm data gives a dust temperature of 7^ = 
23.0 K. This value agrees well with that of Clemens, Yuo & Heyer 
who fit a temperature of = 23.4 K. The model fitted to the IRAS 
data agrees wifli flie SCUBA data at 850 ixm (Visser et al. 2001) 
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freq/GHz 

Figure 3. L675: Upper limits are shown at 2.7 and 4.85 GHz from the Ef- 
felsberg 100 m telescope and the GB6 survey respectively. These data are 
uv matched upper limits. Data at 14.2 to 17.9 GHz from AMI are shown, 
as are points at 850 /im (353 GHz) from Visser et al. (2001; 2003); and at 
100, 60, 25 and 12 from CYH91. A modified greybody spectrum with 
J3 = 1.5 is fitted to the 100 and 60 ^tm data as described in the text and a 
MC spinning dust model from DL98 is fitted to the AMI data points. The 
combined SED is shown as a soUd line. 



assuming the smaller angular size of l'.5, which may be a result of 
the 2 arcmin chop on the SCUBA jiggle map. Using the source size 
of CYH91, this model predicts a flux from thermal dust emission 
at 16 GHz of only 0.22 ^Jy. 

We fit three models to the AMI data points along with the 
upper limits at 2.7 and 4.85 GHz. The first, ignoring the lower 
frequency constraints, is the best fitting power law to the AMI 
channel data alone. The spectral index is not well constrained: 
o^AMi = —0.7 ± 0.7. Including the constraint given by the upper 
limit at 4.85 GHz requires a ^ —1.15. The third model we fit is not 
a power law but the spinning dust spectrum for cold, dark molec- 
ular clouds (Drain & Lazarian 1998; hereinafter DL98). Using a 
source dimension of Ad — 5.5 arcmin from the AMI map, this 
model fits the data with Af(H2) = (2.4 ± 1.3) x 10^5 m'^ (where 
N{K2) = 0.5iV(H)i). Visser et al. (2001) find a peak and averaged 
column density of 9 and 3.4 x 10^^ vcC^ for L675, consistent with 
this value. We note that the extent of L675 in the SCUBA map is 
small (~ 1 '.5), if this is a consequence of the SCUBA chop then the 
coluirm density values may be underestimated. The full spectrum 
for the spinning dust model, along with the best-fit thermal dust 
model, is shown in Fig. 3. 

One further possibility for L675 may be that we are seeing 
an ultracompact Hii region. A spectral index of a = —1.15 would 
be broadly consistent with that expected for such a region. How- 
ever, to have a turn over frequency greater than 20 GHz (assuming 
an electron temperature of ~8000K) would require an emission 
measure in excess of 1.2 x 10^ pccm^^. An emission measure of 
this magnitude would imply a mass for the Hii region of hundreds 
of solar masses. This is in contradiction to Visser et al. (2001) who 
calculate the mass of L675 to be 0.2 M© , making it the lowest mass 
object in our sample. 
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Figure 4. L709: CLEANed AMI combined channel data at 16 GHz are 
shown as contours with the levels as in Fig. 1. NVSS data at 1.4 GHz are 
shown in grey scale and are saturated at 10 mjy peak to emphasise the struc- 
ture in the map. The AMI pointing centre is indicated by a cross and the 
FWHM of the primary beam as a soUd circle. The synthesized beam of AMI 
towards L709, 2'.5 x 2'A, is shown in the bottom left corner as a filled ellipse. 
This map is not primary beam corrected. Example apertures are shown to 
indicate the division between sources "B" and "C" for flux extraction. 



4.2 L709: 

There is no microwave source coincident with the pointing cen- 
tre towards L709, see Fig 4. The field is dominated by a bright 
radio source to the north-east of the pointing centre, denoted "A" 
in Fig. 4. Source "A" lies slightly outside the FWHM of the AMI 
primary beam and has a primary beam corrected continuum flux 
of Sig = 245 ±13 mJy. This is consistent with fluxes from the lit- 
erature of 5o.365 = 374 ±24mJy (Texas), Si ,42 = 270±27mJy 
(Effelsberg 21cm), S2.695 = 310±31mJy (Effelsberg 11cm) and 
'5'4.85 = 330 ± 29 mJy (GB6) indicating a source with spectral in- 
dex a = 0.07 ± 0.02, see Fig. 5. The NVSS flux density for this ob- 
ject is slightly lower, 5'nvss,i.4 = 218 ± 6.6 mJy suggesting that the 
source is extended. A second, and fainter, object in the field is just 
north-west of the pointing centre. Source "B". This object has a flux 
density of 5i6 = 6.0 ±0.5 mJy in the AMI map and is most proba- 
bly a combination of three unresolved point sources, which may be 
found in the NVSS catalogue, all falling within an AMI synthesized 
beam at this position. These sources have a combined flux density 
of 5i 4 = 11.0±0.5mJy, indicating a decrease in flux at 16GHz. 
The slight extension to the south of this object, "C", is not coin- 
cident with any NVSS point sources and may be associated with 
L709. The object is present at the 5(7 level in the combined channel 
map with a flux of Sig = 1.61 ±0.31 mJy. In the consitituent chan- 
nel maps the source is present with flux densities varying between 
1 and 4 a. At this level of significance it is difficult to fit a reli- 
able spectral index, however the flux density appears to be steeply 
falling with increasing frequency. As the source appears point-like 
this would suggest that Source "C" is not associated with L709 but 
is instead a faint steep spectrum extragalactic point source. 



' http : / /www . astro . princeton . edu/~draine/ dust/ 
dust . mwave . spin . html 



4.3 L860: 

At 16 GHz we see a ridge of emission towards L860, see Fig. 6. We 
investigate its spectral properties using analysis Case (1). Although 
it possesses no obvious counterpart in the unsampled CGPS dataset. 
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Figure 5. Radio spectrum of the L709 point source "A". These data points 
are described in the text. 
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Figure 6. L860: CLEANed combined channel AMI data at 16 GHz are 

shown as contours with the levels as in Fig. 1. uv sampled CGPS data at 
1 .4 GHz are shown in grayscale. The AMI pointing centre is indicated by a 
cross and the FWHM of the primary beam as a solid circle. The synthesized 
beam of AMI towards L860, 2!6 x 2M , is shown in the bottom left corner as 
a fiUed ellipse. This map is not primary beam corrected. Example apertures 
for flux extraction are shown to illustrate the division between the regions 
denoted A-E. 



a uv matched image shows the same structure at 1.4 GHz. To the 
north of the field two point-like radio sources may be found (A 
& B), whilst the ridge of extended emission that runs north-south 
across the pointing centre may be divided into three distinct sub- 
regions of emission (C, D & E). The morphology of these sub- 
regions is not well described by a Gaussian model. Although both 
C & D might be considered to be associated with L860 the derived 
flux densities, see Table 2, and their spectral indices indicate that 
there is no excess emission present at microwave frequencies for 
these sources. 



4.4 L917: 

No radio emission can be seen directly towards L917 at 16 GHz, 
although a ridge of emission runs north-south slightly to the west 
of the pointing centre, see Fig. 7. We investigate the spectral prop- 
erties of this ridge primarily using analysis Case (1). This ridge 




freq/GHz 

Figure 8. Above: Radio spectrum of L917 region "A". Below: Radio spec- 
trum of L917 region "C". Data points are uv matched flux densities at 1.4, 
2.7 and 16 GHz from the CGPS, Effelsberg 100 m telescope and AMI, re- 
spectively. See text for details. 



has three separate peaks, which are evident in both the 1.4 GHz uv 
sampled data and the AMI data at 16 GHz. We fit for the flux den- 
sity of each peak separately using the flux extraction method de- 
scribed in Section 4.1. These peaks (A, B & C) all appear to have 
sUghtly more flux at 16 GHz than at 1.4 GHz, see Table 2. Using 
additional data at 2.7 GHz from the Effelsberg telescope uv sam- 
pled under Case (2) we can fill in more of the flux spectrum. The 
peaks A and C show a spectrum consistent with a region of opti- 
cally thin free-free emission, see Fig. 8. In "B", the closest peak 
to the pointing centre we see a large excess at 16 GHz relative to 
1.4 GHz, although it is not clear if this excess is caused by anoma- 
lous emission. At 2.7 GHz the einission has a largely different mor- 
phology making it difficult to constrain a spectrum using the three 
frequencies. A second factor which makes emission from spinning 
dust unlikely in this case is the spectral index across the AMI chan- 
nels, Ofj'^ 2 — 0.54 ±0.36. With a falling spectrum this region is not 
a good candidate for a spinning dust einission, even considering a 
correction for flux loss. It can be seen from the IRAS data at all 
four frequencies that there is a ridge of dust which is spatially as- 
sociated with regions A and B, see Fig. 7. However, in spite of this 
there is no formal IRAS association with L917. The change in mor- 
phology between the radio frequencies could be accounted for by 
this region being composed of a number of small Hii regions with 
different turn-over frequencies into the optically thin regime. 
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Table 2. Flux densities for compact sources found within the AMI fields which have counterparts at 1.4 GHz. 
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4.5 L944: 

The field around L944 contains a number of compact radio sources 
which can be seen in both the AMI data and the corresponding 
NVSS map, see Fig. 9. Five sources near to the pointing centre all 
have non-thermal spectral indices between 1.4 and 16 GHz consis- 
tent with those extrapolated from the WENSS survey at 327 MHz. 
One source is found in the AMI map just north of the pointing cen- 
tre which possesses no counterpart in the NVSS data. This source is 
coincident with the outflow from L944 which is extended towards 



the North and has a flux density of S\(, = 3.5 ± 0.3 mJy in the com- 
bined channel AMI map. 

The spectrum of the this source appears relatively flat across 
the AMI band, aj^ I = —0.07 ± 1 .04, and may be a thermal source 
which remains optically thick to high frequencies turning over be- 
low 16 GHz. It has no counterpart in either the NVSS data at 
1.4 GHz, uv sampled Effelsberg data at 2.7 GHz or the GB6 data 
at 4.85 GHz. It is not possible to constrain the spectral index of the 
possible optically thick component, as the source flux at 16 GHz 
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Figure 7. L917: (a) Combined channel AMI data at 16 GHz are shown as contours with the levels as in Fig. 1. iiv matched CGPS data at 1.4 GHz are shown in 
greyscale. (b) AMI data at 16 GHz are shown as contours with the levels as in Fig. 1. 100 /jm IRAS data are shown in greyscale. The AMI pointing centre is 
indicated by a cross and the IJWHM of the primary beam as a solid circle. These maps are not primary beam corrected. The AMI synthesized beam, 2!8 x 2^ 1, 
is shown as a filled ellipse in the bottom left comer. Example flux extraction apertures are shown to illustrate the divisions between the regions denoted A-C. 
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Figure 9. L944: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. NVSS data at 1.4 GHz are 
shown in greyscale and are saturated at lOmJy peak. The AMI pointing 
centre is indicated by a cross and the FWHM of the primary beam as a solid 
circle. This map is not primary beam corrected. The synthesized beam of 
AMI towards L944, 2'.5 x 2'.3, is shown in the bottom left comer as a fiUed 
ellipse. 



into the cm regime. The second of these seems unlikely: protoplan- 
etary discs are associated with later T-Tauri stars and are also ex- 
tremely smaU in angular diameter. For the flux density measured at 
16 GHz to be accounted for by the tail of the dust SED we would 
require /3 =0.1, consistent with the SED of a protoplanetary disc. 
However, we find that the flux densities from SCUBA and IRAS 
are best fitted by j3 = 0.8 in agreement with Froebrich (2005). As- 
suming a source point-like to the IRAS beam, this would predict a 
flux density from thermal dust of only 0.55 mjy at 16 GHz. 

The possibility of cm-wave emission from a protostar is not 
unlikely. L944 is known to be a Class source, having extended 
and centrally peaked submillimetre continuum emission, which in- 
dicates a spheroidal envelope; a high ratio of submillimetre to bolo- 
metric luminosity (Froebrich 2005) and bipolar molecular outflows 
(Visser et al. 2002). Given these conditions it should also have an 
internal heating source, which would give rise to compact cm- wave 
emission. Previous observations of candidate protostars in the cm 
regime (e.g. Andre et al. 1993; Stamatellos et al. 2007) found fluxes 
which would not be inconsistent with that measured by AMI to- 
wards this source, however we note that the observed AMI flux 
densities would be higher than the average flux density. 



already lies below the noise level for both the Effelsberg data and 
the GB6 survey. 

The short frequency lever arm of the AMI data alone makes 
an estimate of the spectral index from 14.2 to 17.9 GHz imprecise. 
However, including the constraint provided by the lack of any flux 
at 1.4 GHz in the NVSS data constrains the spectral index to be 
aj 4^ < —1.2. Such a spectral index could be consistent with the 
optically thick regime of an ultracompact/hypercompact Hii re- 
gion, the spectral indices of which are normally slightly shallower 
than a = -2 (Franco et al. 2000). 

Two obvious possibilities to explain these observations might 
also be cm- wave emission from a protostar within L944, or alterna- 
tively emission from a protoplanetary disc which may also extend 



4.6 L951&L953: 

There is no radio emission evident towards either L951 or L953, 
see Fig. 10 and Fig. 11. The two fields overlap by a small amount 
- source "D" in the L951 field is also source "A" in the L953 field. 
The flux of this source agrees well between the observations and 
although it appears to have a rising spectrum from 1.4 to 16 GHz 
it has turned over by 4.85 GHz and may be found in the GB6 cat- 
alogue with a flux of 54 gs = 36 ± 5 mJy, giving a spectral index 
of gj = 0.41 ± 0. 14. The ridge of extended emission which runs 
roughly north-south along the eastern edge of the L951 field and 
the western edge of the L953 field may be associated with an en- 
hancement of the local infrared emission and may be seen in IRAS 
data towards the same region, see Fig. 10 (b). 
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Figure 10. L951: (a) CLEANed combined channel AMI data at 16GHz are shown as contovirs with the levels as in Fig. 1. CLEANed uv matched CGPS data 
at 1.4 GHz are shown in greyscale. (b) AMI data at 16 GHz are shown as contours with the levels as in Fig. 1. 100 /im IRAS data are shown in greyscale. The 
AMI pointing centre is indicated by a cross and the IJWHM of the primary beam as a soUd circle. The AMI synthesized beam, 2^5 x 2^2, is shown as a filled 
ellipse in the bottom left comer. 
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Figure 11. L953: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. NVSS data at 1.4 GHz are 
shown in greyscale and are saturated at lOmJy peak. The AMI pointing 
centre is indicated by a cross and the I^WHM of the primary beam as a soUd 
circle. The AMI synthesized beam, 2(5 x I'.'i, is shown as a filled ellipse in 
the bottom left comer. 
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Figure 12. L1014: CLEANed combmed channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1 .4 GHz are shown in greyscale. The AMI pointing centre 
is indicated by a cross and the FWHM of the primary beam as a solid circle. 

The AMI synthesized beam, 2!5 x 2M, is shown as a filled ellipse in the 
bottom left comer. Example flux extraction apertures are shown to illustrate 
the division between sources "A" and "B". 



4.7 L1014: 

The L1014 field contains a number of compact sources none of 
which is completely coincident with the SCUBA position of the 
dark cloud, see Fig. 12. The two closest sources "A" and "B" 
both have falling spectral indices and are therefore not candidates 
for anomalous emission. Like L944, L1014 is a Class source 
(Young et al 2004; Bourke et al. 2005). Radio cm-wave obser- 
vations with the VLA have identified a candidate for the proto- 
star at RA 21h24m07.53s, 8 = +49°59'08.9" (J2000.0) (Shirley 
et al. 2007). Their measurements would imply a flux density of 
5i6 = 145^iJy at 16 GHz, well below the thermal noise in this ob- 
servation. 



4.8 L1021: 

This relatively empty field contains only one radio source which 
can also be identified in the lower frequency data, see Fig. 13. Its 
falling spectral index suggests that it is a non-thermal extragalactic 
point source. 

4.9 L1103: 

This field contains a number of small diameter sources, the bright- 
est of which, "A", we identify as J214307+564033. Sampling the 
corresponding CGPS data towards LI 103 reveals low frequency 
cotmterparts to many of these, see Fig 14. However, the object 
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Figure 13. L1021: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1.4 GHz are shown in grey scale. The AMI pointing centre 
is indicated by a cross and the I^WHM of the primary beam as a solid circle. 
The AMI synthesized beam, 2'.4 x 2(1, is shown as a filled elUpse in the 
bottom left comer. 
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Figure 14. LI 103: CLEANed combined channel AMI data at 16 GHz are 
shown as contovirs with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1.4 GHz are shown in grey scale. The AMI pointing centre 
is indicated by a cross and the I^WHM of the primary beam as a solid circle. 
The AMI synthesized beam, 2(4 x 2(1, is shown as a filled eUipse in the 
bottom left comer. Example flux extraction apertures are shown to illustrate 
the division between those sources labelled C-E. 



found closest to the pointing centre of this field has no lower fre- 
quency counterpart. The maximum flux towards the pointing centre 
from the uv sampled CGPS data at 1 .4 GHz is 1 mjy beam" ^ , com- 
pared to a peak flux density of 4 mJy beam^' at 16 GHz. Moreover, 
an extended structure is obvious in the AMI data running south 
from the pointing centre which is not evident at 1.4 GHz. The total 
integrated flux of this object from the combined AMI channels is 
40. 1 ± 3. 1 mJy. The ridge is present at all AMI frequencies and has 
a rising spectral index of a = —1.6 ±0.3. However the ridge has 
dimensions of approximately 10' x 3', and at such sizes we would 
expect a significant amoimt of flux loss across the AMI band. This 
flux loss would have the effect of steepening a falling spectrum, or 
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Figure 15. Top: SED of L1103. Data points are shown at 1.4 and 2.7 GHz 
from the CGPS and Effelsberg 100 m. These data are uv matched to the 
data at 13.8 to 17.6GHz from AMI which are also shown. As are points at 

850 ^im (353 GHz) from Visser et al. (2001; 2003); and at 100, 60, 25 and 
12^m from Clemens, Yuo & Heyer (1991). A modified greybody spec- 
trum with /3 = 1.5 is fitted to the 850, 100 and 60 ^m data as described 
in the text and a MC spinning dust model from DL98 is fitted to the AMI 
data points, see text for details. The combined SED is shown as a solid 
line. Bottom: SED of J214307+564033. Data points are shown at 327 MHz 
(WENSS, Rengelink et al. 1997), 1.4 GHz (NVSS, Condon et al. 1998; uv 
matched CGPS, this work), 2.7 GHz {uv matched Effelsberg 100 m, this 
work), 4.85 GHz (GB6, Gregory et al. 1996) and 16 GHz (AMI, this work). 
A best fitting spectral index of a = 1 . 1 zb 0.26 is shown as a dashed line. 



flattening a rising spectrum. Since the morphology in the 1.4 GHz 
data is very different we cannot use this as a model to calculate ab- 
solute flux loss. Instead we model the ridge as a multivariate Gaus- 
sian and calculate a flux loss correction for channels 4—8, relative to 
channel 3. This will give us the true slope across the AMI band, if 
not the total power normalization. With these corrections we find a 
corrected spectral index of a = — 2.2± 0.2. A further correction we 
might make is to remove the flux of the point like object "B" from 
the flux density of the ridge as a whole. However, it is unclear from 
a comparison of the AMI and uv sampled CGPS data how large 
the contribution of B to the ridge flux density may be. Although 
the ridge certainly extends to meet B, there is no clear peak at the 
position of the 1.4 GHz source. Consequently, since the difference 
in flux density across the AMI band will be small compared to the 
necessary increase in the error budget we chose not to subtract the 
flux of this object from that of the ridge as a whole. Instead we con- 
sider object "B" as part of the ridge and simply find the integrated 
flux density of the entire ridge. At 1.4 GHz and 2.7 GHz there is no 
apparent elongation towards the pointing centre and so the flux den- 
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sity is heavily dominated by tiie soutiiern end of the ridge. Using 
uv matched data at 1.4 and 2.7 GHz we find a poorly constrained 
spectral index of (X^l = 0.37 ± 0.30. Again we urge caution when 
considering the 2.7 GHz flux density as it is only marginally above 
the noise level of the original Effelsberg data. 

Although, like L675, LI 103 has no formal IRAS association 
it was also observed using deeper IRAS photometry by Clemens, 
Yun & Heyer (1991) and we use their measurements to constrain 
the FIR SED. They found the size of LI 103 to be Q.s = 77 arcmin^, 
more consistent with the area seen by AMI than with SCUBA. 
Again as in the case of L675, this is perhaps not surprising given 
the SCUBA 120" chop. Note that the greybody spectrum which 
best fits the IRAS data points at 60 and 100 flm significantly un- 
derestimates the SCUBA flux at 850/im when the differing angular 
size of the observed source is taken into account. We find a best 
fitting SED with j8 = 0.9 and = 37.2 K, in good agreement with 
fliat of CYH91 who found = 37.6 K with J3 fixed at 1. 

The large extension of the AMI object makes emission from 
either a protoplanetary disc or a single protostar unlikely, although 
the high dust temperature suggests that this is not a prestellar core 
as has been previously supposed. A spinning dust model with the 
dimension of the object fixed at Sis = 30 arcmin-^, as suggested by 
the AMI data would imply aniV(H) = (4.6±0.8) x 10^^ m'^. This 
spectrum is shown in Fig. 15. 

Visser et al. (2001) comment that the distribution of material 
in LI 103 looks as though it may have been swept up by a wind. 
The emission from objects with stellar winds and stars experiencing 
mass loss has been modelled in the past (Calabretta 1991; Casassus 
et al. 2007) by treating the optical depth as a power law function of 
radius. These models, parameterized correctly, can lead to a rising 
flux density spectrum even into the mm-wave regime. We obtain a 
spectral index of OS^gc^™ = — 1.14±0.13, although we note that 
there has been no correction for flux losses in the flux density mea- 
surements used to calculate this index. Indeed, with no basis for a 
model for this source it is not possible to adequately constrain this 
index in order to make any quantitative comments. We note, how- 
ever, that a radial power law model is a vast over-simplification in 
this case as the emission is far from spherically symmetric, and 
also that models of this type which fit both the AMI and SCUBA 
data points greatly over estimate the flux densities measured at 12 
and 25 flm. In addition, these winds are generally associated with 
more evolved protostars. LI 103 is supposedly starless, having no 
observed outflows or IRAS PSC associations. However, the 12 and 
25 jUm flux densities of CYH91 would indicate a hotter dust com- 
ponent to this cloud which might suggest a Class object. 



4.10 Lllll: 

The subject of a separate publication, Lllll possesses a significant 
excess of emission relative to that found at 1.4 GHz in CGPS data 
uv sampled to match the AMI coverage (AMI Consortium: Scaife 
et al. 2009). Using our uv sampling Case (2) we also confirm that 
there is an excess relative to the 2.7 GHz flux density. This map has 
significantly larger resolution (5? 1 x 4'.9) and so the point source to 
the west of LI 1 1 1, denoted "B" in Scaife et al., cannot be clearly re- 
solved from the main object. The integrated flux density of the two 
objects at 2.7 GHz is 15.6±1.0mJy. The spectrum of source "B" 
was measured to be a = 0.43 ± 0. 18, which would give it a flux of 
52.7 = 7-6+} 3 mjy, leaving the object associated with Lllll a flux 
density of S2J = 8.0^} jmjy. However, the noise on the original 
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Figure 16. Lllll: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1 .4 GHz are shown in grey scale. The AMI pointing centre 
is indicated by a cross and the I^WHM of the primary beam as a solid circle. 
The AMI synthesized beam, 2'.4 x 2!l, is shown as a filled ellipse in the 
bottom left corner. Example flux extraction apertures are shown to illustrate 
the division between the sources labelled "A" and "B". 
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Figure 17. SED of Lllll. uv sampled data points are shown at 1.4 and 
2.7 GHz from the CGPS and Effelsberg 100 m. Data at 13.8 to 17.6 GHz 
from AMI are also shown. As are points at 850 ^ra (353 GHz) from Visser 
et al. (2001; 2003); and at 100, 60, 25 and 12 ^tm from the IRAS point 
source rejects catalogue. A power law is fitted to the data at 1.4 and 2.7 GHz 
which is shown as a dashed line. A modified greybody spectrum with J3 = 
1.5 is fitted to the 850 and 100 /im data as described in the text and a MC 
spinning dust model from DL98 is fitted to the AMI data points, see text for 
details. The combined SED across the AMI band is shown as a solid Une. 



Effelsberg data is very similar to this at 8 miy beam" , so the uv 
sampled flux density should be viewed with caution. 

Although there is no IRAS association with Lllll, we fit a 
greybody spectrum to the 100 jim flux from the IRAS RSC Reject 
Catalogue. Using the 100 and 850 ]im flux densities we find a best 
fitting temperature of = 20 K using Q.^ = 6 arcmin^ measured 
from the SCUBA map. The spectral index across the flux loss cor- 
rected AMI data is a\l^ = -2.9 ±0.4, see Scaife et al. 2009 for 
further details. 
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Figure 18. L1166: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1.4 GHz are shown in grey scale. The AMI pointing centre 
is indicated by a cross and the I^WHM of the primary beam as a solid circle. 
The AMI synthesized beam, 2(6 x 2!2, is shown as a filled elUpse in the 
bottom left comer. 
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Figure 20. LI 185: Radio spectrum of Source "A". Data points are at 
365 MHz (Texas Survey, Douglas et al. 1996), 1 .4 GHz (uv matched CGPS, 
this work) and 16 GHz (AMI, this work). A best fitting spectral index of 
a = 1.3 zb 1.1 is shown as a dashed line. 

that derived from data at 365 MHz (Texas Survey) and the CGPS 
at 1.4GHzof a^ |g5 = 1.6± 0.6, see Fig. 20. We therefore suggest 
that this is a steep spectrum extragalactic source. 
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Figure 19. LI 185: CLEANed combined channel AMI data at 16 GHz are 
shown as contours with the levels as in Fig. 1. CLEANed uv matched data 
from the CGPS at 1.4 GHz are shown in grey scale. The AMI pointing centre 
is indicated by a cross and the FWHM of the primary beam as a solid circle. 
The AMI synthesized beam, 2(8 x 2(2, is shown as a filled eUipse in the 
bottom left comer. 



4.11 L1166: 



This object is coincident with a patch of diffuse emission seen 
at 1.4 GHz which has become significantly fainter at 16 GHz, see 
Fig. 18. 



4.12 L1185: 

The L1185 field is relatively empty, see Fig. 19. The compact 
source present to the south of the pointing centre in the uv sampled 
CGPS data has a flux density of 5i.4 = 66.1 ± 1.3mJy at 1.4GHz 
and 5i6 = 2.6± 1 .6 mjy (2cr) in the AMI data. This gives it a poorly 
constrained spectral index of aW = 1.3 ± 1.1, but is consistent with 



4.13 L1246: 

LI 246 is a reasonably well studied dark cloud. It is known to con- 
tain a protostar (SMMl; Launhardt et al. 1997; Visser et al. 2001), 
IRAS 23228+6320, at 23h25m05.6s +63°36'34".l (J2000.0). The 
position of the IRAS association seems to be dominated by the 25 
and 60 ^im emission, with the 12 and 100 pLm maps showing emis- 
sion further to the east. The AMI map at 16 GHz shows a band of 
emission which traces that found at 12 ^m, see Fig.21. The north- 
ern most lobe of which is also coincident with an 8 /im point source 
seen in the MSX' image of the same region, see Fig. 21. This emis- 
sion has no counterpart in the lower frequency uv sampled CGPS 
data at 1.4 GHz. The known protostar L1246 SMMl can also be 
seen in the 8 ^um data just to the west of the pointing centre. Two 
of the four PAH mid-infrared (MIR) features at 7.7 and 8.6 /im 
lie within MSX Band A. The unidentified MSX point source co- 
incident with the AMI emission in the L1246 field is only present 
within this band. 

We note that although Visser et al. (2001) take the size of this 
cloud to be small, using the dimensions of Parker (1991) of 2!8 x 
\'A, the cloud appears larger in the DSS plates with dimensions 
much closer to those measured by Clemens & Barvainis (1988) of 
lOf 1 X 2!2. It is unsurprising therefore that in a later paper Clemens 
et al. (1991) measure a much larger flux density towards this object 
than that listed for the IRAS association of Parker (1988). Given 
the large angular extent of the region in the co-added IRAS images 
of Clemens et al. (1991) it seems likely that these flux densities will 
come from an area covering the emission seen in the AMI data. 

Previous greybody fits to L1246's SED have suggested = 
19 K (Launhardt et al. 1997, Visser et al. 2001) using data at A > 
100 /im. The temperatures found using IRAS data alone are higher. 
From the IRAS PSC the temperatures found using 100 to 60 /im, 60 
to 25 /im and 25 to 12 /im are 31, 66 and 101 K respectively. From 
the deep IRAS photometry of CYH91 using a much larger aperture 
the same temperatures are 31, 69 and 222 K. These fits are sensitive 



' http : / /irsa . ipac . caltech . edu/applicat ions/MSX/MSX/ 
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Figure 21. L1246: AMI 16 GHz combined cliaimel data are shown as contours at 3, 4, 5, 6, 9, 12 and 15 a. Greyscale is (a) uv matched COPS data at 1.4 GHz, 
(b) 8 lira MSX (Band A) data, (c) 100 /jm IRAS data, and (d) 12 /im IRAS data. In plot (a) the position of the protostar SMMl is indicated and in plots (c) 
and (d) the position of the associated IRAS point source is shown as a square. In all plots the AMI pointing centre is shown as a cross and the primary beam as 
a solid circle.. The AMI synthesized beam, 2'.5 x 2'.2, is shown as a filled ellipse at the bottom left of each image. The intensity in the MSX and IRIS images 
is saturated at both ends of the scale in order to make the structure in this region visible. 



to the angular size assumed for the object, the discrepancy at 12 jUm 
suggests a significant change in morphology at that wavelength. 

The morphology of the AMI data at 16 GHz appears to trace 
the 12 /im emission and to a lesser extent the 100 flm emission also. 
Taking a 5 arcmin square region around the emission seen in the 
AMI data to the east of the pointing centre there is a slightly higher 
spatial correlation between the emission seen at 12 nm {r = 0.70) 
with that seen at 16 GHz than there is at 100 /im (r = 0.61). There 
is no source coincident with the AMI emission at either 1.4 or 
2.7 GHz. The AMI data itself towards this source is poor due to 
interference and the combined channel map is heavily weighted 
towards Chaimel 6 (16.4 GHz) making a determination of the spec- 
trum using AMI data alone unreUable. 

5 DISCUSSION AND CONCLUSIONS 

Of the fourteen clouds observed here two are clear candidates for 
possessing anomalous dust emission at 16GHz: L675 and LlUl. 



A further three are possible candidates: L944, L1103 and L1246. 
Nine of the fourteen showed no emission inconsistent with that seen 
at lower radio frequencies. 

We divide the five clouds where there is an excess of emis- 
sion into two groups on the basis of three criteria: (1) Extent of the 
emission seen at 16 GHz i.e. non-point like sources; (2) Coincident 
with the sub-mm position of the cloud to within a 2' radius; (3) 
Likelihood of alternative explanations for the 16 GHz emission. 

We summarize our results in Table 3. Coluirm [5] of Table 3 
summarizes which dark clouds have a possible spinning dust asso- 
ciation, defined as being a microwave counterpart within a 2 arcmin 
radius which shows an excess at 16 GHz relative to lower frequency 
data. Columns [6-8] then divide the certainty of these detections 
based on the three criteria stated in the previous paragraph. 

On the basis of these three criteria we see that all are satis- 
fied by L675 and LI ill which display extended emission coinci- 
dent with the pointing centre. The extended nature of this emis- 
sion makes the possibility of the emission being cm-wave radia- 
tion from either a protostar or proto-planetary disc unlikely. L944 
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shows compact emission which is coincident with the pointing cen- 
tre. The compact nature of this emission means that cm-wave radi- 
ation from a protostar cannot be ruled out in this case. L1246 shows 
emission just within a radius of 2 arcmin from the SCUBA position 
of this cloud which has no lower frequency counterpart. However, 
the compact nature of this emission suggests that it may again be 
cm-wave emission from a protostar. 

LI 103 shows extended emission which is coincident with the 
pointing centre. However, although the microwave spectrum ap- 
pears to be consistent with emission from spinning dust this source 
has a high flux density at 850 /xm (Visser et al. 2001) relative to 
IRAS measurements at 12-100 ^im (CYH91). This raises questions 
over the suitability of a standard greybody fit to the thermal dust 
emission from this cloud and, in the absence of further observations 
between 100 and 300 GHz, it is not possible to rule out a flattened 
tail to the thermal dust spectrum which may account for the excess 
seen at 16 GHz. 

In their theoretical models for the emission from rotating 
grains DL98 considered a number of gas-grain, plasma, and ra- 
diative contributions to the excitation and damping of spinning 
dust grains. For their several models of grain environment it was 
found that collisions with ions were the dominant excitation mech- 
anism in all but two scenarios: reflection nebulae and photodisso- 
ciation regions, where photoelectric effects dominated due to the 
intense radiation field. For the environments relevant to dark neb- 
ulae (molecular cloud and cold neutral medium enviromnents) this 
suggests that objects which are embedded in or exposed to strongly 
ionized environments are the most likely candidates for observing 
spinning dust emission. The presence of an anomalous component 
in the spectra of the objects presented here appears to be correlated 
with the average free-free emission within a 1 degree area centred 
on each pointing centre, see Fig 22. We assess this correlation fol- 
lowing the method of Franzen et al. (2009), where we express the 
covariance matrix of our data points, C, as 



S+Ni 



Sii S22 
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(5) 
(6) 



5i 1 and 522 are the variances of the 2.7 and 16 GHz data points 
respectively andSi2(=52i) is the covariance of the 2.7 and 16 GHz 
data points. It follows that the correlation co-efScient, denoted R to 
distinguish it from the standard Pearson co-efflcient, is given by 



R- 



S12 



VSuS22 



(7) 



Following Franzen et al. (2009) we used a Gaussian likelihood 
and the nested sampling algorithm implemented in the MultiNest 
code (Feroz & Hobson 2008; Feroz, Hobson & Bridges 2008) to 
obtain R. For the averaged 2.7 GHz flux density we find a corre- 
lation coefficient of i? = 0.81 ± 0. 15. For the 2.7 GHz flux density 
at the position of the cloud with the 16 GHz flux density at the 
same position we find a weaker correlation, R = 0.66 ±0.20. We 
repeated these evaluations setting the noise covariance matrix to 
zero and found values oiR = 0.70 ± 0. 18 and = 0.65 ± 0.20 for 
the averaged flux density at 2.7 GHz and the 2.7 GHz flux density 
at the position of the cloud, respectively. We can see that the noise 
in the first instance leads to an decrease in correlation unless it is 
accounted for. 

Of the fourteen clouds observed, those with the highest back- 



ground levels at 2.7 GHz are L675, LI 111, L860 and LI 103; all 
of which have a background level above lOOmK. Indeed the cor- 
relation is dominated by the high surface brightness regions, sug- 
gesting a radio surface brightness threshold between those regions 
which exhibit spinning dust emission and those which do not. 

There appears to be no correlation between 16 GHz specific 
intensity and 100 ^um emission. Following the same correlation 
analysis a described in Equ. 7 we obtain a value of = 0.02 ± 0.32 
(R = —0.01 ±0.32 without noise). The dust emissivity at 16 GHz 
relative to the IRIS 100 nm data is tabulated in Table 3. We also 
list the predicted specific intensities for these objects for a relative 
emissivity of 10 /iK(MJy/sr)^' for reference. The relative emis- 
sivity of our two strongest detections is similar to that found for 
LDN1622 (Casassus et al. 2006) and G159.618.5 (Watson et al. 
2005). This suggests that the non-detections in our sample are not 
an effect of our sensitivity limits. If the emissivity of anomalous 
microwave emission relative to that seen at 100 /im was at a simi- 
lar level in all dark clouds, for example > 10^iK(MJy/sr)~', then 
we would have detected emission from every cloud in our sam- 
ple at > 3<Jrms- The weighted average dust emissivity relative to 
the 100 /xm maps for our sample is 6.6±5.0 /iK(MJy/sr)^'. This 
value agrees with the all-sky WMAP value for cool dust clouds at 
high galactic latitudes (Davies et al. 2006), given as 1 Jy at 30 GHz 
per 6000 Jy at 100 /im. 

Of our five possibly anomalous objects two are known to be 
Class objects: L944 and L1246. High sensitivity maps in CO to- 
wards our two most anomalous objects do not yet exist and they are 
assumed starless. However the starless/protostellar divide has come 
under scrutiny recently following the discovery of Very Low Lumi- 
nosity Objects (VeLLOs) with Spitzer. VeLLOs are faint infrared 
point sources with protostellar colours towards objects previously 
classified as starless. Young et al. (2004) discovered a VeLLO in 
L1014 (L1014-IRS), part of the Visser et al. sample, coincident 
with the peak of dust continuum emission that had been previously 
classified as starless by Parker (1988). Although here no anomalous 
microwave emission is seen towards L1014 it is suggestive that the 
peak of the anomalous emission towards L944 is more coincident 
with the centre of the red-shifted CO 2^1 outflow lobe, mapped 
by Visser et al., and not the compact SCUBA emission from the 
Class protostar. This spatial correspondence immediately raises 
the question of whether the outflow itself or possibly the cloud tur- 
bulence could be energising the spinning dust emission. Regardless 
of this possible connection we also point out that the misidentifi- 
cation of radio emission from protostars is a potential source of 
confusion when searching for emission from spinning dust. 

There seems to be no obvious correlation between the AMI 
flux densities and the IRAS or SCUBA flux densities, however the 
number of data points is low. It is also true that there is no corre- 
lation between the SCUBA flux densities at 850 /xm and those of 
IRAS at 100 /xm. This is perhaps not particularly surprising in light 
of the mismatch in measured angular scales, and would suggest 
that the morphology of these objects is larger at 850 /tm than the 
SCUBA chop allows for. Obviously this is not a problem when in- 
vestigating compact cores as in Visser et al. (2001; 2002), however 
it provides little assistance when looking at the correlation between 
cm- wave amd mm-wave emission. 

In conclusion, we have observed a sample of fourteen compact 
Lynds dark nebulae. We have found a significant excess towards 
two of the fourteen and an indication of anomalous behaviour in 
three further clouds. We suggest that the excess we see is due to ro- 
tational emission from very small grains and that this emission may 
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Specific Intensity at 1 6GHz [mJy beam ] 

Figure 22. Correlation of mean radio brightness temperature with microwave specific intensity towards observation pointing centres. Data are specific intensity 
at 16 GHz towards the SCUBA position of the dark cloud on the x-axis, and mean brightness temperature within a 1° diameter circle centred on the position 
of the dark cloud on the y-axis. The error bars on this mean show the standard deviation of the brightaess temperature distribution within this region. 



liable 3. 16 GHz dust emissivities relative to 100 /im emission. Column [1] Source name. Column [2] Predicted intensity based on an emissivity of 
10 ;UK (MJy/sr)~^ , Column [3] measured intensity towards phase centre, Column [4] Measured 16 GHz dust emissivity relative to IRIS 100 ^m data, Col- 
aam [5] Possible identification of a spinning dust counterpart as described in Section 5, Colvmm [6-8] Certainty criteria as described in Section 5. 



Name 


(mJy beam"') 


5act 

(mJy beam"') 


e 

/iK(MJy/sr)"l 


Possible 
Identification 


[1] 


Criterion: 
[2] 


[3] 


L675 


3.6 


5.9±0.3 


16.4±1.8 


y 


y 


y 


n 


L709 


2.8 


-0.8±0.3 


<0 


n 








L860 


4.2 


3.0±0.3 


7.2±1.0 


n 








L917 


7.9 


-0.7±0.3 


<0 


n 








L944 


2.7 


1.7±0.3 


6.2±1.3 


y 


n 


y 


n 


L951 


1.9 


-0.9±0.3 


<0 


n 








L953 


2.4 


1.1±0.6 


4.5±2.5 


n 








L1014 


2.5 


0.0±0.4 


0.0 


n 








L1021 


2.8 


1.0±0.4 


3.5±1.4 


n 








L1103 


4.7 


1.8±0.5 


3.9±1.2 


y 


y 


y 


y 


Lllll 


3.1 


8.7±0.5 


28.5±3.3 


y 


y 


y 


n 


LI 166 


1.5 


0.3±0.6 


2.0±4.0 


n 








LI 185 


4.2 


-0.4±1.1 


<0 


n 








L1246 


2.1 


0.2±0.4 


1.0±2.0 


y 


n 


y 


n 



be correlated with a high level of background emission at lower ra- 
dio frequencies. 
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